Critical currents at the Bragg glass to vortex glass transition 
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We present simulations of the transport properties of superconductors at the transition from the 
Bragg glass (BG) to the vortex glass (VG) phase. We study the frustrated anisotropic 3D XY 
model with point disorder, which has been shown to have a first order transition as a function of 
the intensity of disorder. We add an external current to the model and we obtain current-voltage 
curves as a function of disorder at a low temperature. We find that the in-plane critical current has 
a steep increase at the BG-VG transition, while the c-axis critical current has a discontinous jump 
down, this later result in agreement with the first-order character of the transition. 
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The study of the vortex phase diagram of anisotropic 
superconductors with point disorder has been of interest 
since the discovery of high Tc superconductivity. It is now 
clear that at low temperatures and low magnetic fields 
there is a "Bragg glass''_phase (BG) with an elastically 
distorted vortex lattice |1| . This vortex lattice undergoes 
a first order melting transition to a vortex liquid (VL) 
when increasing the temperature T. At higher magnetic 
fields and low T there is a disordered vortex state, the 
"vortex glass" (VG) A disorder driven transi- 

tion from the BG to the VG has been proposed I"] to oc- 
cur when increasing the magnetic field H or the disorder. 
Experimental observations with increasing H at low T, 
that have been attributed to the BG-VG transition, are: 
the destruction of the Bragg peaks in neutron scattering 
0, a dip in the differential resistance 0, the onset of 
the "second magnetization peak" or a jump in the 

Josephson plasma resonance |lOj . Numerical simulations 
in particle- like models with random pinning have found a 
transition from an ordered lattice to a disordered lattice 
[n, JIJ when increasing particle density (i.e., the mag- 
netic field). Recently, Monte Carlo simulations in the 
frustrated 3D XY model with disorder have stablished 
that the BG-VG transition is of first-order type 13, nMj. 

Several of the experimental evidences of the BG-VG 
transition are directly or indirectly related to transport 
properties. For example, the most common determina- 
tion of the BG-VG line is through the onset of a "sec- 
ond magnetization peak" 0] , which is attributed to an 
steep increase of the in-plane critical currents. From the 
point of view of transport, the BG has zero linear resis- 
tivity (piin = 0) The VG phase was originally pro- 
posed to have zero resistivity 1|. However, studies of the 
XY gauge glass model with finite screening {k < 00) have 
found that VG is unstable in d = 3 and the disor- 
dered phase is a frozen vortex liquid that has a very small 
finite resistivity. Only for the k — > 00 disordered XY and 
gauge glass models the VG phase is stable in = 3 at low 
T, in which case pun — 0^. Therefore, it is important to 



understand how the the non-linear current-voltage (IV) 
curves and their critical currents change across a disor- 
der driven /irsi-orrfer transition. We present here numer- 
ical calculations of both the in-plane and the c-axis IV 
curves across the disorder driven BG-VG transition in 
the K = 00 frustrated XY model used in [HQ (which 
is the only model where the first-order nature of the tran- 
sition has been clearly determined so far). We find that 
the c-axis critical current has a discontinuous jump at the 
BG-VG transition, which could be tested experimentally. 
The hamiltonian of the frustrated 3D XY model is: 

where (r) = 9{r + p) - 9{r) - A^(r) and A^(r) ^ 

^ /J^^^"*" A-dl and r (nx,ny,nz) defines a lattice site 
in a cubic grid. To minimize the pinning of the numerical 
grid, we have chosen V{6^{r)) ~ —Rq — Ri cos(6'^(r)) — 
R2 cos(20„(r)) and adjusted the coefficients Ro,Ri and 
i?2 as in We model uncorrelated random point pin- 
ning in the xy plane |l3l ll^ with J^f^ = J{1 + per,fi), 
/i = X, y, where er,^ are independent random variables 
with (er,^) = and (e^^^J = I. In the presence of an 
external magnetic field H along the z direction, we have 
Ax{r) - Ax{r + y) + Ay{r + x) - Ay{r) = 27r/, and we 
take / — Ha^/^o = ^. We model the current between 
two grid points with the RSJ model [l^ : 

I^{r,t)=Sf,{r,t) + N^{r,t)+ 7]^ (r , t) 

with the superconducting current ^^(r, t) — x ae^(r} ' 

normal current N^{r,t) — 2e/? ^^at''^ with i?^ — Rat for 
/i = X, y and R^ — TRab for p, — z, and the thermal 
noise fluctuations have correlations (7/^(r, t)?/^/ (r', t')) = 
{2kBT / R^j)5jx^jx'5r.r'5{t — t'). We take periodic bound- 
ary conditions in all directions with a fluctuating twist 
such that A^{Y,t) = Aj^(r) — Q;^(i). This allows to 

obtain the voltage in each direction as — 55—^^ 



2 



[l7| . We consider current conservation in each node: 
^/x(r) — Ifii''^ — fJ-) = 0: and we fix the total current in 
each direction, consistently with the periodic boundary 
condition, by |l7l |: 



rcxt 



2eR^~dr 



1 



LxLyLz 



5][5^(r)+7y^(r,i)] 



We integrate with a second order Runge-Kutta 
stochastic algorithm with At = 0.1 h'^ / [Ae^ RabJ) and we 
average over an interval of 10^ — 10^ At. We consider sys- 
tem sizes 48 X 48 X with Lz — 12 — 32 and anisotropy 
r = 40. Most of the results shown are for = 32 and for 
a single disorder realization. We have considered other 
realizations of disorder for some values of T and p and 
have found similar behavior. 
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FIG. 1: (a) Intensity of the Bragg peak S(G) and c-axis re- 
sistivity Pa = Vz/h vs. r for p = 0.01 and size 48 x 48 x 32. 
T is normalized by J/fcs. (b) 5(G) vs. p for T = OAJ/ks- 
A: 48 X 48 X 12; ■: 48 x 48 x 32. Left inset: intensity plot of 
S(k) for p = 0.07. Right inset S(k) for p = 0.14. 



We first study the model at equilibrium (in the absence 
of external currents). At low disorder, p = 0.01, we find 
a transition upon increasing T from a vortex lattice to 
a vortex liquid. We have calculated the structure factor 
^(k) = (l/iV/2)(|Er,,.".(r±,^)e^''-'""n, where is 
the vorticity in the xy planes and N = L^LyLz- In 



Fig. 1(a) we plot the intensity of the structure factor, 
S{G), at one of the first reciprocal lattice vectors G, 
showing the melting of the vortex lattice at Tc « 0.28. 
The resistivity is calculated by driving the system with 
a very small current and obtaining voltage. We show 
in Fig. 1(a) the c-axis resistivity pc vs. T which drops 
sharply at T^, we find similar behavior in pab (not shown). 
At a low temperature T = OAJ/ks, we vary disorder and 
study the behavior of the lattice order parameter, 5(G), 
which is shown in Fig. 1(b) for two system sizes. We find 
that >S'(G) drops at pc ~ 0.11, indicating a transition 
from the BG to the VG. Typical structure factors S'(k) 
for p < Pc and p > pc are shown in the insets. Repeating 
this analysis for other T, we have found that Pc{T) is 
approximately independent of T for low temperatures. 
Similar phase boundary Pc{T) for the BG phase were 
reported in For very low temperatures, both 

and agree that there is a first-order transition from 
the BG phase to the VG phase in this model. (At higher 
T there is disagreement on the existence of a "vortex 
slush" phase). 

From now on, we focus on the disorder driven transi- 
tion at low T, where [l^ll^ agree on the phase diagram. 
We fix a low temperature T = O.lJ/kB and we vary the 
disorder strength p. In Fig. 2(a) we show the calculated 
in-plane IV curves for different p. They are obtained 
by applying a uniform current along the x direction, 
jcxt _ I^^Sfj^x, and calculating the average voltage along 
X, Vab = (dax/dt). In Fig. 2(b) we show the c-axis IV 
curves. In this case we apply a uniform current along the 
z direction (parallel to the magnetic field), 1°^*" = Iz5^i,z: 
and we calculate the voltage Vz — (daz/dt). Before do- 
ing further analysis it is important to discuss the effect 
of a small driving current in the BG-VG transition. We 
find that currents much smaller than the critical cur- 
rents have a negligible effect on pc- In Fig. 3 we show 
5(G) vs. p when the driving current in a given direction 
is of the order of the corresponding critical current. In 
the case of a z-direction current Iz — 0.010 ~ there 
is a very small effect with a critical Pc slightly smaller 
than the 7 = case. In the case of an in-plane current 
lab = 0.002 ~ 7°^, we find that there is a first drop of 
5(G) to a finite value at a p < Pc- This corresponds 
to the situation when the applied current lab equals the 
critical current for that p. Then a second drop of 5(G) 
to the very small values of the VG phase occurs near the 
critical pc but also at a slighlty smaller value. In both 
cases, one can conclude that finite currents of the order 
of the critical current have a small effect on the value of 
Pc, lowering it in about Apc ~ —0.005. 

Let us now discuss the in-plane transport at the tran- 
sition. In Fig. 2 (a) we show with open symbols the IV 
curves for p < pc and with filled symbols the IV curves 
for p > Pc (where Pc = 0.11 is the zero current critical dis- 
order). We observe that in both phases the IV curves are 
strongly nonlinear and that they tend to zero resistivity 
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FIG. 2: (a) In-plane lab-Vab curves for different p. (b) c-axis 
Iz-Vz curves. (/ normalized by Jo = 2eJ /h and V by Ratio)- 



FIG. 4: (a) In-plane voltage Vat vs. p. (b) c-axis voltage Vc 
vs. p. (Voltages normalized by Ratio, size is 48 x 48 x 32.) 
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FIG. 3: ^(G) vs. p for finite currents, and 48 x 48 x 32 size. 



for low currents {V/I — > 0). In order to observe if there 
is a change when crossing the BG-VG line, we fix a value 
of the current lab, slightly above the critical currents, 
and we calculate the voltage Vab as a function of disorder 
strength p. This is shown in Fig. 4(a). We see that at 
Pc there is a jump down in Vab- This jump is observed 
for different values of lab- The location of the jump does 
not depend strongly on lab- The experiment of [7| finds 



a clear decrease of the differential resistance, measured 
for a current above the critical current, when increasing 
magnetic field at low temperatures in YBaCuO samples. 
This is in agreement with our result of Fig. 4(a) for the 
behavior of the in-plane voltage at the BG-VG transi- 
tion. One can interpret this jump down in voltage as a 
consequence of a jump up in the critical current. It is 
difficult to define a critical current, since at T > there 
is always a non-zero voltage for any current. Experimen- 
tally, the "critical current", Ic, is obtained from the IV 
curves as V{Ic) = Vth, with the threshold voltage Vth 
small enough. If 1^ is almost independent of Vth, then it 
corresponds to a crossover current separating two regimes 
in the IV curve. In our case, Ic separates the ohmic flux 
flow regime from the low current nonlinear regime which 
has p\in{I ^ 0) = for the two phases. (In the case of the 
VG with flnite screening, k < oo, 1^ would separate two 
ohmic regimes). To obtain the ah critical current, I^^, 
we choose different values of Vth slightly above our noise 
value of voltage. One of the Vth used is shown in Fig. 2(a) 
with a dashed line. The obtained critical currents are 
plotted in Fig. 5 (a) as a function of p (we see that differ- 
ent Vth give smilar values of if' , within the error bars). 
We observe in the plot that there is a steep increase of 
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/° at Pc- In the magnetization experiments in YBaCuO 
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and we calculate 14, as we show in Fig. 4(b). For very low 
Iz we find 14 « for all p, consistent with both phases 
being superconducting. For intermediate values of Iz we 
find that Vz ~ for p < pc with a sharp jump to finite 
voltage Vz at Pc- For values of Iz above the critical cur- 
rents, there is a finite voltage Vz in both phases, with a 
clear discontinuous jump up at pc- We obtain the critical 
currents I^ with a voltage criterion Vth, which is shown 
in Fig. 2(b) with a dashed line. It is obvious in this case 
that the values obtained for I^ are insesitive to the choice 
of Vth ■ The c-axis critical currents are plotted as a func- 
tion of disorder strength p in Fig. 5(b). We observe in 
the plot that I^ has a clear discontinuous jump down at 
Pc- This result is consistent with the Josephson plasma 
resonance measurements of , which are sensitive to 
the c-axis Josepshon coupling, where a sharp decrease 
of the plasma resonance frequency was observed at the 
transition. Recently, the c-axis Josephson critical cur- 
rent has been measured experimentally in BiSrCaCuO 
samples varying magnetic field at the melting transition 
line, the BG-VL line A jump down of the critical 

current was found, although at a field slightly lower than 
the equilibrium case. It will be interesting if a similar 
experiment could be carried out across the BG-VG line, 
since our results imply that c-axis transport would give 
a better evidence (compared to the a6-plane transport) 
of the first-order nature of this transition. 



FIG. 5: (a) In-plane critical current J"'' vs. p, for threshold 
voltages Vth- (b) c-axis critical current vs. p. 

m the BG-VG transition line has been associated to the 
onset of the second magnetization peak, which has been 
related to an steep increase of the in-plane critical cur- 
rents. Our result in Fig 5(a) is in good agreement with 
this interpretation, particularly because the value of F we 
use corresponds to YBaCuO. The more anisotropic BiS- 
rCaCuO samples show experimentally a much stronger 
increase in the second magnetization peak 0. In this 
very anisotropic case the increase in critical current has 
also been interpreted as a decoupling transition . 

The 3D XY model is very adequate to the descrip- 
tion of the c-axis transport, since the Josephson coupling 
between planes is treated exactly. The c-axis current- 
voltage charactestics, Vz vs. Iz, are shown in Fig. 2(b) 
for different values of p at T = O.lJ/fcs. The IV curves 
for p < Pc are plotted with open symbols and the IV 
curves for p > Pc are plotted with filled symbols. We ob- 
serve that in all cases the IV curves show a well-defined 
critical current where the voltage drops steeply. The IV 
curves can be clearly separated in two different sets of 
curves, one set for p < Pc with high critical currents and 
another set for p > pc with low critical currents. This 
can also be observed if for a fixed current Iz we vary p 
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